引言

  1. 定义:特征值全是正数的实对称矩阵为正定矩阵(positive definite matrix)
  2. 类似的,若实对称矩阵的特征值均非负,则为半正定矩阵(positive semidefinite matrix)。

可能用到的概念

  • 主子式

    • 定义:在nn阶行列式中任选kk行,再取相应的kk列,将行列交汇处元素组成新的矩阵行列式,称为nn阶行列式的一个kk阶主子式。
  • 顺序主子式(the k-th leading principal minor)

    • 定义:在nn阶行列式中由第1,...,k1,...,k行和第1,...,k1,...,k列所确定的主子式称为kk阶顺序主子式。直观上看就是矩阵中左上方的子矩阵。

实对称矩阵AA正定的充要条件

注意,这里的所有进行判别的矩阵都是实对称矩阵。

以下条件都是判别实对称矩阵是否正定的充要条件。

  1. AA的所有特征值λi\lambda_i均为正。
  2. xTAx>0x^TAx > 0对所有非零向量xx都成立。
  3. AA的所有顺序主子式都是正的。
  4. AA的所有主元(无行交换)都是正的。
  5. 存在列满秩矩阵RR,使得A=RTRA=R^TR.
  6. AA的所有主子式都是正的。

证明

(1) => (2)

对实对称矩阵AA,存在正交阵QQ,使得A=QΛQRA=Q\Lambda Q^R

因此,对任意非零向量xx

xTAx=xTQΛQTx=yTΛy=λ1y12+...+λnyn2>0x^TAx = x^TQ\Lambda Q^Tx = y^T\Lambda y=\lambda_1y_1^2 + ... + \lambda_ny_n^2 > 0

其中y=QTx=(y1,...,yn)0y = Q^Tx = (y_1,...,y_n) \ne 0

(2) => (1)

因为AA为实对称矩阵,因此一定可以对角化A=QλQTA = Q\lambda Q^T

xTAx>0xTQΛQTx>0\because x^TAx > 0 \quad \therefore x^TQ\Lambda Q^Tx > 0

y=QTxy = Q^Tx,因此yΛyT>0λiyi2=0y\Lambda y^T > 0\Rightarrow \sum \lambda_iy_i^2 = 0

y0s.tλiyi2=0\because \forall y \ne 0 \quad s.t \quad \sum \lambda_iy_i^2 = 0

λi=0\lambda_i = 0

(2) => (3)

由(2) => (1) => det(A)=λ1...λn>0det(A) = \lambda_1 ...\lambda_n > 0

xTAx=(xkT0)(Ak)(xk0)=xkTAxk>0x^TAx = \begin{pmatrix}x_k^T & 0\end{pmatrix}\begin{pmatrix}A_k &* \\\\ * & *\end{pmatrix} \begin{pmatrix}x_k\\\\0 \end{pmatrix} = x^T_k A x_k > 0

det(xkT)det(Ak)det(kk)=det(Ak)>0\therefore det (x_k^T)det(A_k)det(k_k) = det(A_k) > 0

(3) => (4)

顺序主子式与主元有直接关系,第kk个主元:

dk=det(Ak)detAk(k1)>0d_k = \frac{det(A_k)}{detA_k(k-1)} > 0

其中AkA_k是第kk个顺序主子矩阵(the k-th leading principal submatrix)。

(4) => (2)

由对称矩阵的Gauss消元法得LDU分解:A=LDLTA = LDL^T,其中对角阵的对角元为AA的主元:

D=diag(d1,...,dn)D = diag (d_1,...,d_n)

则对任意非零向量:

xTAx=xTLDLTx=yTDy=d1y12+...+dnyn2>0x^TAx = x^TLDL^Tx = y^TDy = d_1y_1^2 + ...+ d_ny_n^2> 0

(2) => (5)

A=LDLT=LDDLT=(DLT)T(DLT)=RTRA = LDL^T = L\sqrt {D} \sqrt {D} L^T = ( \sqrt {D} L^T)^T ( \sqrt {D} L^T) = R^TR

(5) => (2)

A=RTRA = R^TR,则对任意非零向量xx:

xTAx=xTRTRx=(Rx)T(Rx)=Rx2>0x^TAx = x^TR^TRx = (Rx )^T(Rx) = ||Rx||^2 > 0

(6) => (2)

(6) => (3) => (2)

(2) => (6)

kk阶主子矩阵Ai1,...,ikA_{i1,...,ik},任取x=(x1,...,xn)0x = (x_1,...,x_n) \ne 0,使其除xi1,...,xikx_{i1},...,x_{ik}的其余分量全为0,则

xTAx=(xi1,...,xik)Ai1,...ik(xi1...xik)>0x^TAx = (x_{i1},...,x_{ik})A_{i1,...ik}\begin{pmatrix}x_{i1}\\\\... \\\\ x_{i_k}\end{pmatrix} > 0

det(Ai1,...ik)>0\therefore det (A_{i1,...ik}) > 0

判断

因此,我们可以用以上充要条件来判断矩阵是否正定。常用的判断方法有:

  1. 看顺序主子式是否都大于0
  2. Gauss消元后主元是否都大于0
  3. 看特征值是否都大于0
  4. 找任意一个向量计算xTAxx^TAx是否大于0
  5. A=LDLTA = LDL^T分解,看R=DLTR = \sqrt {D} L^T是否满秩

正定矩阵的常见性质

  1. 正定矩阵一定可逆。(行列式大于0)

  2. A,BA,B为正定矩阵,则A+BA+B也为正定矩阵。

  3. AA为正定矩阵,则存在矩阵C,使得A=C2A = C^2

    • 证明:AA为正定矩阵,则存在正交矩阵QQ,使得:

    • A=QΛQT=(QΛQT)(QΛQT)=C2A = Q\Lambda Q^T = (Q\sqrt{\Lambda}Q^T) (Q\sqrt{\Lambda}Q^T) = C^2

  4. AA为正定矩阵,则矩阵A2A^2A1A^{-1}也正定。

  5. AA为正定矩阵,矩阵CC可逆,则B=CTACB = C^TAC也正定。

半正定矩阵

充要条件

  1. AA的所有特征值λi\lambda_i均非负。
  2. xTAx0x^TAx \ge 0对所有向量XX成立。
  3. 存在矩阵RR,使得A=RTRA =R^TRRR可能不是可逆阵)。
  4. AA所有主子式均非负。(注意,不是所有顺序主子式

results matching ""

    No results matching ""